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Abstract. In this paper we investigate the overcooling 
problem and propose some possible solutions. Wc show 
that the overcooling problem is generic to the hierarchical 
picture of structure formation, as long as the cooling is 
actually possible in small halos at high redshift. Solutions 
to this problem are likely to be associated with some feed- 
back mechanism, and probably imply the existence of a 
warm IGM containing most of the cosmological baryons. 
We concentrate on some possible solutions, mainly pho- 
toionization and bulk-heating of the IGM. We show pho- 
toionization can act as a significant feedback mechanism 
but is not strong enough to solve the entire overcooling 
problem. We therefore assume that the IGM is maintained 
hot by some energy injection provided by supernova and 
galaxy formation is then limited by this feedback mecha- 
nism. Such a self regulated scheme allows us to compute 
the thermal history of the IGM. In the absence of any 
photo-ionization, we find the existence of bifurcations in 
the thermal history. However, these bifurcations are sup- 
pressed when significant photo-ionization occur in addi- 
tion to the energy injection. The temperature of the IGM 
is then maintained at T ~ few 10 5 K. We find that for 
realistic fraction of the energy produced by supernovae 
being re-injected into the IGM, this scenario can consis- 
tently reproduce the present amount of stars. A T~ CDM 
model with T ~ 0.25 can reproduce properly the observed 
HI gas at high redshift. A simple prescription for star for- 
mation also allows to reproduce the star formation rate 
at high redshift as inferred from recent data, while the 
entropy of the IGM at z ~ 2 — 5 is of the order of what 
seems necessary to explain the observed properties of X- 
ray clusters. We conclude that the warm self-regulated 
IGM picture provides an interesting alternative to stan- 
dard semi-analytical approach, which may elucidate the 
behavior of baryons in structures from small galaxies to 
clusters. 
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1. Introduction 

The global picture of the formation of structures on all 
scales is one of the widest questions in modern cosmology. 
Even the validity of the gravitational instability paradigm 
has not yet received complete confirmation and the door 
remains open to alternative approaches, like the topolog- 
ical defect models (see for instance Durrer et al. 1999). 
However, the gravitational instability picture is the most 
compelling one, primarily since it can be investigated to 
an unprecedented level of detail. The small perturbations 
necessary to seed structure formation are believed to orig- 
inate in the very early universe. The obvious non-linear 
character of structures like clusters and galaxies has lead 
to the intensive use of numerical simulations which are 
now achieving remarkable performances. In this context , 
the role of baryons has also been included in recent years 
by incorporating the hydrodynamics of gas in the simula- 
tions. However, analytic modeling is also essential to our 
understanding of the physics, even if some aspects require 
the validation by numerical simulation. The clearest ex- 
ample of this is the mass function as inferred by Press 
and Schechter: it is only since numerical simulations have 
confirmed its validity, at least as an efficient fit to the ac- 
tual mass function arising in numerical simulations, that 
it has been widely used for cosmological applications. Cos- 
mological constraints that can be inferred from properties 
of clusters is certainly the field which has the most bene- 
fited of this approach, although it is not yet entirely clear 
whether we understand the physics of the gas in clusters. 
For instance, the relation between luminosity and tem- 
perature is not what one would expect in a simple scal- 
ing model as proposed by Kaiser (1986): the origin of the 
problem probably lies in the fact that the formation pro- 
cess of the X-ray core is not well understood (Blanchard & 
Silk 1990) possibly requiring non-gravitational processes 



2 



S.Prunet 1 ' 2 & A. Blanchard 3 ' 4 : Warm photo-ionized IGM: a Clue for galaxy and cluster formation history? 



(Evrard & Henry 1991; Kaiser 1991). However, the impres- 
sive agreement of the Press and Schcchtcr formula with the 
mass function on every scale, for which some theoretical 
reasons exist, has lead to handle the question of struc- 
ture formation on galactic and sub-galactic scales up to 
the question of the first structures that formed (Tcgmark 
et al. 1997). On galactic scales it is clear that dissipative 
processes play a key role, as they are necessary to allow 
star formation. It has been proposed that the cooling cri- 
terion is an essential requirement for galaxy formation: 
in order to form a galaxy it seems actually unavailable 
that the baryons can dissipate their thermal energy and 
reach a density high-enough that star formation can ac- 
tually start, although the details of this process are still 
far from being understood: even star formation in our own 
Galaxy or in neighboring galaxies is not well understood. 
The fact that first stars should form in the absence of 
any dust makes the situation more problematic. The in- 
clusion of the baryonic component in PS halos was pio- 
neered by White and Rees (1978). They tried to infer the 
shape of the luminosity function of galaxies, and noticed 
that it was expected to be steeper than the mass function 
at the faint end, leading to a luminosity function far too 
steep compared to observations. This was also noticed by 
Peacock and Heavens (1990) in the context of the peak 
formalism approach to the mass function. The reality of 
this problem has become obvious, in more recent analysis, 
in which several authors have tried to infer the shape of 
the luminosity function in the context of the CDM the- 
ory (White & Frenk 1991; Lacey & Silk 1991). In their 
approach, White and Frenk (1991) developed a model for 
the formation and evolution of galaxies in which a lot of 
energy was injected into the gas which has felt in the dark 
matter potential wells, something which was achieved at 
the price of a metal production level which appeared un- 
realistically large. Since, similar models have been build 
in order to achieve a better agreement with observations 
(Kauffmann et al. 1993; Cole et al. 1994). In the con- 
text of the Cold Dark Matter picture in an Einstein de 
Sitter universe Blanchard et al. (1990; 1992a), BVM here- 
after, investigated the cosmological history of the cooling 
of the gas. They concluded that the simple cooling argu- 
ment leads to the conclusion that most of baryonic gas was 
likely to have been cooled at some time of its history, lead- 
ing to the so called overcooling problem: the amount of gas 
which would be expected today to be in cold and probably 
dense state would represent more than 80% of the cosmo- 
logical baryonic content of the universe, one order of mag- 
nitude larger than the baryonic content of galaxies. This 
problem was also noticed independently by Cole (1991). 
The overcooling problem is a great challenge to numer- 
ical simulations as spatial resolution of 10 3 is necessary 
and a mass resolution of 10 6 (first significant cooling oc- 
curs in objects which have a virial temperature of 10 4 K at 
Z ~ 10 corresponding to mass of 10 8 M Q with typical size 
of 1 kpc). However, the high resolution numerical simula- 



tions of Navarro and Steinmetz (1997) were able to see this 
overcooling problem, with an amount of cooled gas which 
agrees quite well with a simple argument based on Press 
and Schcchtcr estimate. This gives enough confidence in 
the analytical approaches for further investigations of the 
problem, and its possible solutions. Indeed, any reliable 
scenario for galaxy formation should solve this problem 
in a self-consistent way. BVM proposed that the IGM is 
heated by the feedback of first structures to temperature 
high enough (typically a few 10 5 K) to prevent further 
collapse in most potential wells. In this picture most of 
the IGM stays homogeneous. This warm IGM picture was 
shown to potentially offer an elegant solution to the ori- 
gin of the shape of the luminosity function, leading to an 
inverse hierarchical picture, in which large galaxies form 
at high redshift while at least some small galaxies form at 
low redshift. One of the objectives of the present work is to 
investigate this idea in a more quantitative way. Further- 
more, several new observations provide information rele- 
vant to galaxy formation: the CFRS redshift survey has 
illustrated the apparent increase of the star formation rate 
up to redshift of the order of one, while the HDF seems 
to reveal a decrease of the star formation rate at higher 
z. Dust obscuration may however occur at large redshift, 
masking a significant fraction of star formation. Further- 
more, the first field galaxies at redshift as high as 4 have 
been discovered (Steidel 1996) and the neutral gas content 
in dense objects has been evaluated up to redshift 4, this 
gas potentially being the precursors of disks. It is there- 
fore interesting to see if such observations can be accom- 
modated by this model. In the first section, we examine 
the overcooling problem in various models and its sensi- 
tivity to the various input parameters. We then examine 
the possible role of photo-ionization and show that photo- 
ionization implies a significantly smaller fraction of cooled 
gas, but is still unable to solve the overcooling problem. 
We then investigate in some detail the original scenario 
of BVM in which the IGM is heated by the feedback of 
galaxy formation, resulting in a self regulated evolution 
of the temperature of the IGM. This model will be called 
the warm IGM picture hereafter. The basic conclusion of 
our study is that the warm IGM picture possesses few pa- 
rameters but can easily describe several global properties 
of the cosmological baryonic history that can be inferred 
from the observations. It may therefore represents an in- 
teresting alternative to more traditional semi-analytical 
scenarios. 

2. The Overcooling problem 

The physics of the formation of galaxies is still rather un- 
known. However, it is almost trivial to say that the physics 
of baryons should play an essential role. For instance, the 
origin of the characteristic size, mass and luminosity of a 
bright galaxy has been first looked for in the physics of 
the linear evolution of fluctuations before recombination 
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or after the recombination (see for instance Gamov 1948). 
However, the Jeans masses after and before recombina- 
tion were recognized as being very different from the typ- 
ical galaxy mass, while the damping leads to masses too 
large to be identified with galaxies (Silk 1967). One key 
criterion which has been proposed to define the scale asso- 
ciated with typical galaxies (say an L* galaxy) is based on 
the cooling criteria: during the gravitational collapse of an 
object the kinetic energy of the gas is turned into thermal 
energy, which is then lost by some mechanism in order for 
the gas to contract further up to densities where star for- 
mation has any reasonable chance to start (of course, this 
is a minimal requirement, as the mechanism of star for- 
mation may require many more conditions!). It has been 
suggested that galaxies are the largest objects able to 
cool in a time scale shorter than the typical age of the 
universe (Binney 1997; Rees & Ostriker 1977; Silk 1977). 
These kinds of considerations were already investigated by 
Hoyle as early as 1953 (Hoyle, 1953). 

It is in this context that the overcooling problem has 
been outlined, which is likely to occur in any hierarchical 
model of galaxy formation. The nature of the overcooling 
has been explicitly stated for the first time by Blanchard 
et al. (1990; 1992a) and Cole (1991). In the simplest pic- 
ture, one assumes that the gas which cools is turned into 
stars, so the total amount of cooled gas should represent 
the present day amount of stars, this is quite reasonable 
as the amount of cold gas (HI) at low rcdshift is modest. 
In the context of hierarchical models, it is expected that 
at some redshift a large fraction of the matter of the uni- 
verse settles into small potential wells with circular veloc- 
ity between 20 and 200 km/s, the virial temperature of the 
gas being then at temperature higher than 10 4 K and at 
rather high density. Under these conditions the radiative 
cooling is very efficient and therefore the cooling time is 
extremely short. One then naturally expects a large frac- 
tion of the cosmological gas to be cooled by now. This is in 
contradiction with a standard result in Cosmology: stan- 
dard nucleosynthesis value for the baryonic density pa- 
rameter of the universe is of the order of 0.08/i^ 2 (Buries 
et al. 1999), while the present density in stars There are 
several considerations which are important to take into 
account while dealing with the overcooling problem: us- 
ing the cooling criterion as a criterion for star formation 
reflects our poor knowledge of the realistic physical con- 
ditions required for star formation. However, the cooling 
criterion seems clearly to be a minimal requirement, and it 
is likely to be valid anyway in order to evaluate the total 
amount of cold gas, counting stars in this category. Be- 
cause the present amount of cold dense gas actually seen 
in the universe is rather small, it is reasonable to imagine 
that cold gas has been actually turned into stars. Some 
alternatives are possible though: there might be more gas 
than what we see (Pfenniger & Combes 1994), or most of 
the cold gas may have been turn into stars which are not 
seen. In both cases, this represents a way to generate can- 
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didates for baryonic dark matter. In such a case, the over- 
cooling question would be intimately linked to the solution 
of the dark matter problem! This could at least point out 
to a possible origin of the machos in galactic halo. Machos, 
if they exist, may provide evidence that a substantial frac- 
tion of the baryonic density of the universe is dark, at least 
if the situation in our Galaxy can be taken as representa- 
tive of the general situation (see Alcock et al. 1998). One 
could also wonder if the apparent overcooling problem is 
not due to the fact that nucleosynthesis estimates are erro- 
neous. However, this can not provide a satisfying answer. 
The reason comes from the x-ray properties of clusters : in 
cold dark matter picture the baryonic content is expected 
to be the sum of the stars in galaxies and of the x-ray 
gas, the latter representing at least 80%/i~ 3//2 of the total 
baryonic content, meaning that the fraction of gas which 
is in the "cool" phase (stars in this case) is a small fraction 
of the total baryonic content. More precisely, the typical 
ratio of the baryonic mass seen in clusters to the mass in 
stars is of the order of 5-6 h~ 3 / 2 (White et al. 1993), so 
assuming that this ratio is universal leads to an estimate 
of the baryonic content of the universe, mainly in the form 
of intergalactic gas : 

ttbaryon^cluster « (l + 5-6ft- 3/2 )ft, ~ 0.03-0.06/l- 3/2 (l) 

interestingly close to primordial nucleosynthesis value ~ 
0.02ft, -3 / 2 for a low Deuterium abundance. Incidentally, 
this has interesting consequences: this value conflicts with 
nucleosynthesis values corresponding to high Deuterium 
abundance (Songaila et al. 1997) and does not allow a 
significant contribution of machos to the cosmic baryon 
budget. The agreement between the inferred value of fib 
and nucleosynthesis reinforce the so-called baryons crisis 
(White et al. 1993), which is a problem specific to fl = 1 
models. Finally, the ratio of mass in visible stars to gas 
mass in clusters is consistent with the ratio F t , the ratio 
between 0* and Clbbn in the universe: 

F* - 0.05 - 0.20 (2) 

This essentially means that only a small fraction of the 
primordial baryons were turned into visible stars and that 
there is little room for an other significant dark baryonic 
component. 

2.1. Framework 

There has been a lot of progress in the investigation of 
cosmic scenarios for structure formation since it has been 
realized, thanks to numerical simulations (Efstathiou et 
al. 1988; Cen, Gncdin & Ostriker 1993; Navarro, Frcnk 
& White 1996), that the Press and Schcchter (1974) for- 
mula, (PS hereafter), for the mass function is a good work- 
ing approximation. Our understanding of the behavior of 
the baryonic content has also been greatly improved by 
numerical simulations in which the fate of the cosmologi- 
cal gas is followed. Despite possible complex behavior, it 
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Fig. 1. In this figure we have estimated the integrated 
fraction of gas F c (using Eq ||) in two T— CDM pictures 
with r = 0.5 (thin lines) and T = 0.25 (hick lines) with 
different values of the bias parameter with rj — 1 . As one 
can see the history of the cooling depends on the param- 
eters of the model, but the present day total amount of 
baryons that have cooled is relatively independent of the 
model, and reaches a large value, close around 80%. 

seems well established that the gas is heated during the 
collapse either by shocks or by adiabatic compression, un- 
til it reaches the virial temperature at which it can settle 
in nearly hydrostatic equilibrium: 

r„~5x 10 5 M^ 3 (l + z) K (3) 

However, hydrostatic equilibrium is possible only as long 
as the gas does not loose a substantial fraction of its in- 
ternal energy, i.e. during a period of time smaller than 
the typical cooling time. The PS mass function changes 
over a time scale comparable to the age of the universe, 
so if the typical cooling time is longer than the age of the 
universe, the gas will be merged in an other structure be- 
fore it can cool. On the contrary, if the gas can cool on a 
time scale much shorter than the age of the universe, it 
is unlikely that gravitational shaking will re-heat the gas 
before it has cooled down. As soon as the gas can cool it 
will contract further. The gas will be re-heated during the 
contraction phase by adiabatic compression, however as 
the cooling time becomes shorter with increasing density 
there is a runaway catastrophe, which can apparently be 
stopped only when the gas has dissipated its thermal en- 
ergy, probably being then rotationally supported or having 
been turned into stars. 

The instantaneous fraction of gas able to cool can be 
computed as a simple integral: 

1 r™ 2 

f c {z) = — N(m)mw(m)dm (4) 
P Jmi 

where w(m) represents the fraction of the baryonic ma- 
terial present in the structure m which is able to cool or 



Fig. 2. On this figure we have computed the instanta- 
neous fraction of gas f c in two T— CDM pictures with 
r = 0.5 and T — 0.25 with different values of the baryonic 
content of the universe with r/ = 1 and b = 1.5. 

has already cooled at an earlier epoch. BVM showed that 
w{m) can be evaluated by assuming that all the bary- 
onic content is in a homogeneous phase inside the halo 
and estimating its cooling time: if this cooling time t c (m) 
is shorter than some characteristic time t\, the gas will 
cool before undergoing any substantial re-heating. Such a 
characteristic time scale will always be a fraction of the 
Hubble time at the considered redshift: 

h = i]tH (5) 

with rj smaller than 1, and that can be taken to be in- 
dependent of the mass. This may not be strictly valid: 
the typical time scale may depend on the power spectrum 
and on the amplitude of the initial fluctuations associated. 
However, it is doubtful that the effect is important. Ac- 
cordingly to BVM, we assume that the cooling fraction is 
either 100% or 0% i.e. w(m) can take only two values (0 
or 1) and cooling is efficient only for masses in some range 
mi — ni2- The instantaneous cooling fraction can then be 
written as: 

f c (z) = - mN(m)dm (6) 

P J mi 

The next issue is to estimate the total fraction of cooled 
gas integrated over redshift. As long as the integrated frac- 
tion remain small one can writes: 

F c (z)* f +X '£&dt (7) 

where rjtjj is a typical time scale for haloes to rearrange 
themselves. When the integrated fraction starts to become 
large this formula cannot be used any more: it becomes 
likely that most of the baryonic content of the halos in the 
cooling region has already cooled in a previous generation 
of halos. In principle, one has to follow the detailed history 
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at all epochs of the precursors of these haloes, which could 
be done for instance by using a merging history tree. How- 
ever, BVM showed that the result of this computation can 
be evaluated more directly by taking into account the past 
history of the halos: large masses at low redshift are essen- 
tially built from structures which were inside the cooling 
region at some higher redshift. Therefore, they should al- 
ready contain only cooled gas. This leads to: 

1 f +oc 

F c {z) = - I N(m)mdm (8) 
P J mi 

avoiding to have to handle the merging tree history. Nu- 
merical simulations indicate that this is a good working 
approximation (Navarro & Steinmetz 1997). In a realis- 
tic scenario however F c (z) has to remain small and the 
approximation ^ can be used. 

Let us now examine the typical amplitude one gets for 
the amount of cooling at various redshifts and the sensi- 
tivity to the various parameters. We have evaluated the 
integrated fraction F c (Eq. ||) for two CDM-likc models 
with shape parameter fj T = 0.5 and T = 0.25 (hence- 
forth referred to as T-CDM models) with various values 
of the normalization. This is presented in figure [y. Here 
and throughout the article we adopted the cooling func- 
tion with zero metallicity as given by Sutherland & Dopita 
1993. This figure illustrates the overcooling in a dramatic 
way: at redshift z = more than 80% of the primor- 
dial baryons should have been cooled. This is almost in- 
dependent on the value of T and on the the normalization 
erg = 1/6, although the redshift history is sensitive to the 
specific model. 

Because the amplitude of the overcooling problem 
is nearly a factor of ten and because the final value is 
rather insensitive to the detailed values of the parameters 
entering the calculation, it is not expected to be solved 
by adjusting one of the various parameters entering 
the problem. Examination of the instantaneous fraction 
(figures H reinforces the idea that the problem is real: 
there are epochs at which the instantaneous amount of 
gas able to cool is in the range 20% - 50% (the maximum 
value of f c gives a robust lower bound for the integrated 
value F c ). This is already larger than what is implied by 
the observed amount of stars: the estimated present-day 
value is F* ~ 10%, Varying the baryonic content of 
the universe or the value of T of the power spectrum 
or the normalization of the models does not change 
significantly the maximum value of f c . This illustrates 
the significance of the overcooling problem. Another 
interesting consequence is that the amount of cold gas 
available at high redshift is high enough that most of the 
models can accommodate a main galaxy formation epoch 
as early as z > 5: in fact in the context of the standard 

1 r, the shape parameter, is denned as the apparent magni- 
tude of Q,h in the CDM transfer function, see e.g. Peacock & 
Dodds 1994 



cold dark matter picture (with b ~ 2), all the observed 
stars could have formed at z > 5 ! The Canada-France 
Redshift Survey (CFRS) has provided evidence that a 
significant amount of star formation has occurred since 
z = 1 (Lilly et al. 1995; 1996) and the Hubble Deep 
Field (HDF) may provide evidence that the bulk of star 
formation occurs at low redshift, possibly as low as z — 1 
(Madau et al. 1996; 1997). A qualitative comparison of 
the instantaneous fraction f c {z) and of the star formation 
rate as inferred from observations is instructive: from 
/ c (z), one expects the bulk of star formation to occur 
at rather low redshift (zt ~ 1 — 4). In other words, 
although the amount of gas available for star formation 
is much too large compared to the SFR, the shapes of 
the redshift distributions of both quantities are not very 
different (although a very late galaxy formation epoch, 
like z < 1, would not appear natural). This is consistent 
with the idea that galaxy formation has occurred in a 
self-regulated way by some mechanism, although this 
mechanism is not explicitly stated in most semi-analytic 
models. As argued by BVM, because such a mechanism 
should essentially limit the amount of star formation, 
it is reasonable to think that the same mechanism is 
responsible for the solution of the overcooling problem 
and might explain the shape of the luminosity function of 
galaxies. The scenario we investigate in the next section 
is an example of such a scheme. 



3. Re-heated IGM and structure formation 

Although there is little doubt that re-heating is a fun- 
damental problem during galaxy formation, it is much 
harder to identify which astrophysical sources could be 
responsible for it, and to specify the detailed physics of 
the feedback mechanism. 

Following earlier consideration by Larson (1974), Dekel 
and Silk (1986) argued that star formation was self limited 
in small potentials by supernova heating: in this model 
gas is heated and expelled from galaxies. The same mech- 
anism was advocated by Cole (1991). This effect is ex- 
pected to be significant for dwarf galaxies with masses 
M & 1O 7 M (see Mac Low & Ferrara, 1999). These au- 
thors also showed that if gas cannot be expelled from 
larger galaxies, metals can still be expelled efficiently for 
masses up to M ~ 10 8 AfQ). Tegmark et al. (1993) have 
reexamined such a model and concluded that the IGM 
could be re-ionized by supernova-driven winds, solving the 
Gunn-Peterson test, but probably producing a substantial 
amount of metals at redshift z = 5 or larger. 

A feedback mechanism has been incorporated into 
most recent models of galaxy formation, although it is gen- 
erally implemented by assuming that galaxy formation is 
inefficient within small halos, while the detailed physics is 
not explicitly taken into account. This approach is the one 
used in several so-called semi-analytic models (eg Kauff- 
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mann et al. 1993; Cole et al. 1994). In our opinion, this 
procedure might hide one of the most critical aspects: it 
is far from obvious that the overcooling problem has been 
solved then, even if the global modeling seems to match 
observations, because the physics of the feedback mech- 
anism is not explicitly treated. There are basically two 
different kinds of feedback mechanisms. The first family 
is the one used by White and Frenk (1991), and is local 
in nature. The second advocated by BVM is global. The 
suppression of cooling by photo-ionization enters this last 
category. Certainly one of the central issues of galaxy for- 
mation is to understand whether the feedback mechanism 
which actually solves the overcooling problem is local or 
global. 

In the re-heated picture, the physical state of the IGM 
might be quite inhomogeneous and complex. In the sce- 
nario which we will consider we simplify the general pic- 
ture by assuming that there are essentially two phases: one 
which is the condensate phase, corresponding to the gas 
which has been able to cool and can eventually be turned 
into stars, and the IGM phase, which is supposed to be 
essentially homogeneous. At high redshift, as soon as the 
IGM is ionized, its cooling time is shorter than the Hub- 
ble time and it lies well inside the cooling phase (unless 
its temperature is very high, but then this would produce 
unacceptably large Compton distortion in the CMB spec- 
trum, given the constraints provided by the COBE data, 
see Fixsen et al. 1997) and some heating process is neces- 
sary to maintain it in this hot phase. 

3.1. Structure formation suppression mechanisms 

In order to solve the overcooling problem one needs to 
satisfy two constraints: the integrated density of stars 
should not be larger than what is observed in present-day 
visible stars and the amount of cooled gas at high redshift 
should not be larger than what has been estimated 
from Damped Lyman alpha systems. Although there are 
ways to escape these constraints, in the following we 
will investigate which scenarios naturally satisfy these 
constraints. 

In order to reduce the amount of cooling that might 
happen during the cosmological history of baryons, there 
are two interesting mechanisms that one can think of. The 
first mechanism has been advocated by BVM: if the IGM 
is hot with a typical temperature Tigm then structure for- 
mation will be suppressed on scales for which the virial 
temperature is smaller than this: 

Tv < ?IGM 

This is certainly a robust criterion, but it is possible that 
the actual suppression is much more efficient: if the IGM is 
at some finite temperature its contraction during gravita- 
tional collapse could be adiabatic and its temperature can 
rise. In a strict adiabatic collapse, one might expect that 



the actual suppression occurs on mass scales for which the 
temperature satisfies the following criterion: 

T v < A 2 / 3 T IGM 

It is interesting that for Tigm ~ 10 5 K, this leads to the 
suppression of structure formation for halos with V c ~ 100 
km/s. However the existence of a cut-off for galaxy for- 
mation at a single circular velocity is not very attractive: 
observations do reveal the existence of small galaxies and 
the luminosity function does not reveal any specific feature 
around M ~ 16. Actually, there might be an increase in 
the number of faint galaxies at this luminosity rather than 
a decrease as one would expect (ESO Slice Project, Zucca 
et al. 1998). Another possible mechanism for suppressing 
galaxy formation is through photo-ionization. Not only 
does photo-ionization heat the gas to a temperature of the 
order of Tigm ~ 10 4 K, but it can also modify the cool- 
ing function. The implications are discussed in the next 
section. 

3.2. Effect of photo-ionization 

The existence of quasars at high redshift is enough to 
ensure that photo-ionization is playing a role in the his- 
tory of the cosmological baryons. Whether the UV flux at 
high redshift is high enough to ensure the low level of HI 
optical thickness observed towards high redshift quasars 
is still a matter of debate (Giallongo et al. 1996; Cooke 
et al. 1997; Devriendt et al. 1998). In any case, photo- 
ionization by QSO's can easily heat the gas at a temper- 
ature of the order of 10 4 K at redshift as high as 5. In ad- 
dition photo-ionization suppresses collisional line cooling, 
by suppressing the existence of neutral atoms (Efstathiou 
1992). It is important to examine whether photoioniza- 
tion can severely limit the collapse of baryons within po- 
tentials as originally proposed by Efstathiou (1992). This 
has been examined in some details by mean of numer- 
ical simulations. Quinn et al. (1996) demonstrated that 
photoionization has essentially no effect on the collapse 
of halos for which the virial temperature is greater than 
10 4 K. We can therefore compute the amount of instanta- 
neous cooling in a totally photo-ionized medium by using 
the cooling function of a fully ionized gas, assuming that 
only bremsstrahlung cooling is important (see for instance 
Thoul & Weinberg 1996). The resulting f c is presented in 
figure |j| By comparing to figure pL one can see that the 
effect of photo-ionization is to substantially suppress the 
amount of cooling occurring at low redshift. The cumu- 
lative fraction of cooled gas is then reduced by a factor 
of the order of 2 by redshift 0. This implies that photo- 
ionization substantially changes the cooling process of the 
gas at low redshift, but can apparently not solve the over- 
cooling itself. This is because the bulk of the overcooling 
occurs at a redshift when the density of the gas is high 
enough that photo-ionization cannot prevent most of the 
cooling. Collin-Souffrin (1991) has also examined the effect 
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Fig. 3. Same quantity as in figure |2j the instantaneous 
fraction of gas f c in two T— CDM pictures with r = 0.5 
(thin lines) and T — 0.25 (thick lines) with different values 
of the baryonic content of the universe with b = 1.5 in the 
photo-ionized IGM. 



Fig. 4. Same quantity as in figure [l] the integrated frac- 
tion of gas F c in two CDM pictures with T = 0.5 (thin 
lines) and T = 0.25 (thick lines) with different values of the 
baryonic content of the universe with r\ = 1 and b — 1.5 
in the photo- ionized IGM. 



of photo-heating by the x-ray background, but this is not 
an efficient mechanism for raising the temperature much 
above few times 10 4 K. Our study indicates that photoion- 
ization does not solve the overcooling problem, i.e. that 
a large fraction of baryons collapse in small potentials at 
redshift below 10. This picture is rather different from 
what numerical simulations seems to indicate: when pho- 
toionization is taken into account the emerging picture is 
that most of the IGM lies in moderately overdense region 
(1 + A ~ 0.1 — 10), and that no significant collapse occurs 
in small potentials. However, in order to see the collapse 
of baryons in the photoionized regime, great care must 
be taken to spatial and mass resolution (Weinberg et al, 
1996). Let us examine this : present day high resolution 
numerical simulations (see for instance Machacek et al., 
1999) achieve typically 256 3 cells with a box of the order 
lOh -1 Mpc corresponding to a spatial scale of 40/i _1 kpc 
(comoving). At redshift z = 3, this is the typical virial 
radius of a halo with a virial temperature of T ~ 2.10 5 K, 
significantly higher than the temperature of photoionized 
gas. This artificial cut-off in the mass implies these numer- 
ical simulations cannot describe properly the formation of 
halos with virial temperature in the range 10 4 — 10 5 K, in 
which most of the cooling take place in the redshift range 
3 — 10. This suggests therefore that these numerical simu- 
lations are still significantly limited by resolution, and for 
this reason do not see the overcooling problem. This is a 
serious problem for the association of Lya clouds to mod- 
erately overdense clouds might be an artifact due to in- 
sufficient resolution: the IGM might well be actually more 
clumpy than what is found in numerical simulations. 

4. Galaxy formation in the warm IGM picture 



4--1. Heating processes 

As photo-ionization seems not to be sufficient to solve the 
overcooling problem, one has to advocate other heating- 
processes. Star formation is the source of a significant 
amount of energy liberation. Not only massive stars would 
radiate UV photons, but they would also end up in super- 
novae which could inject a significant amount of energy in 
the IGM either through direct mechanical energy either 
via cosmic rays. The first possibility has been examined 
by Tegmark et al. (1993). From theoretical considerations, 
it is possible that a significant of energy produced by a 
supernova is transferred to cosmic rays (Malkov & Voelk 
1995). This is also suggested by the observed amount of 
cosmic rays in the galaxy Drury et al. (1989). Prantzos 
& Casse (1994) have discussed possible evidence that the 
flux of cosmic rays was higher at the beginning of galaxy 
formation, in connection with the star formation process 
itself, in order to explain the abundance of boron (see also 
Nath & Biermann 1993 and Ginzburg & Ozernoi 1965). 
Both mechanisms can easily provide heating sources for 
the IGM. Furthermore, cosmic rays may propagate rather 
easily to large distances, and therefore heat the IGM in 
a rather uniform way. It is therefore well possible that 
a substantial fraction of the energy produced by super- 
novae is transferred to the IGM. In the following we will 
therefore assume that a fraction of the energy produced 
by supernovae is transferred to the IGM. The energy from 
supernovae transferred to the IGM will therefore be writ- 
ten in a parametric way: 

Uigm\\+ — £f>*EsN (9) 

where p* ~ \jtu{z) is the star formation rate per unit 
mass, Esn is the energy produced by supernovae per unit 
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mass of stars formed for a standard IMF and e is the 
efficiency of the transfer mechanism. Assuming that the 
cooled gas is instantaneously transformed into stars , the 
heating source of the IGM then is: 



PbUlGAlfc 



(10) 



Assuming a uniform temperature, Tigm, for the IGM con- 
taining all the baryons (i.e. neglecting the few baryons 
which are turned into stars) , the equilibrium temperature 
in the self-regulated case is given by the following equa- 
tion: 



I f" 12 

(Pb/m p ) 2 A(T) = pbUiGM- N{m)mdm 

P J mi 



(11) 



In this equation mi — mi is the range of masses in which 
cooling is able to occur. The lower temperature will be 
determined by the mass of the smallest objects able to 
form in a warm IGM. As we have argued, baryons will 
not collapse in potentials which end up with a virial tem- 
perature smaller than the temperature of the IGM. On the 
other hand, as long as the IGM is in the cooling area, the 
baryons will easily follow the collapse of the dark mat- 
ter potentials provided that the cooling time of the gas 
is shorter than the dynamical time. This means that the 
mass mi corresponds to the mass for which the virial tem- 
perature is equal to the temperature of the IGM. In the 
previous equation, mi therefore is a function of T, and the 
system is fully determined and hence self-regulated: if over- 
cooling occurs a large amount of stars will form, causing 
a large energy injection, thereby raising the temperature 
of the IGM, suppressing galaxy formation by inhibiting 
further collapse. In practice, the system will evolve to- 
ward an equilibrium situation and the temperature can 
be computed. Such a solution may however not necessar- 
ily exist if the cooling is too efficient. This, for instance is 
the case at high redshift, when the number of heat sources 
is limited and the IGM is dense, so that cooling is very 
efficient. 



4-2. Supernovae Feedback 

We will try in the following to get some idea about the 
feed-back parameter e, and about the subsequent produc- 
tion of metals in the IGM. It is known that the progenitors 
of the SNII are stars of mass greater than 8Mq (Rcnzini 
et al. 1993). What we need to compute is the mechanical 
(either thermal or kinetic) energy released per unit mass 
of stars. This will provide E$n, the energy that the SNIIs 
might release into the IGM. We thus need to assume an 
Initial Mass Function (IMF) for star formation, which we 
will take to be the Salpeter IMF: 



6.0 



5.5 



5.0 



4.5 



4.0 



5 10 15 

redshift 

Fig. 5. An example of the equilibrium temperatures of the 
IGM in the self-regulated picture for two specific models 
: the parameter e in two T— CDM pictures with T = 0.5 
(thin lines) and T = 0.25 (thick lines). The baryonic con- 
tent of the universe is set to 0.1, the normalization of the 
power spectrum is b = 1.5 and the feedback parameter is 
set to e = 0.5. 

where x = 1.35 and ty(m) is the differential number of 
stars of mass m. We can then compute the ratio p in mass 
of the SNIIs progenitors to the total mass of stars: 



P 



J s ° m^{m)dm 



Jo^ m\l/(m) dm 



0.21 



(13) 



for a Salpeter IMF. In the same way we compute the av- 
erage mass of the SNIIs progenitors (Msnii) — 30M Q . 
Taking the average mechanical energy released by an SNII 
explosion to be (Esnii) = 10 51 ergs, we can compute €q: 



E, 



SN 



p(Esnii) 
(Msnii) 



3.510 15 ergs/g 



(14) 



*(m) = Am- (1+x) 



(12) 



Only a fraction e of this energy will actually be available 
for heating the IGM. However, as we may consider a flatter 
IMF (x = 1.0) or a bimodal one (Elbaz et al. 1995), it is 
possible in principle for this parameter to be greater than 
one. For this reason, we will consider a range of models 
going from e = 0.125 to e = 2. 

4-3. The warm IGM case 

In this section we will consider the case where only ther- 
mal energy is injected into the IGM. As we already pointed 
out, in our scenario the only free parameter, for a given 
IMF, is e, the efficiency with which energy from super- 
novae is transferred to the IGM. The temperature of the 
IGM is then determined for each specific cosmological 
model. The general equilibrium solutions are presented in 
figure | for some models. Due to the various "peaks" in 
the cooling function several solutions exist for the value 
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Fig. 6. The temperature of the IGM in the self-regulated 
picture for various values of the parameter e in two 
T-CDM pictures with T = 0.5 (thin lines) and T = 0.25 
(thick lines). The baryonic content of the universe is set 
to 0.05 and the normalization of the power spectrum is 
b = 1.5. 

of equilibrium temperature. Some of these solutions are 
unstable and are therefore rejected. One can see that no 
continuous solution exists over the whole range in red- 
shift. In practice, the temperature has to "jump" from 
one equilibrium solution to another, a process known as a 
bifurcation. In such case, the exact behavior of the tem- 
perature can be followed only by numerical integration 
of the non-equilibrium equations and is likely to depend 
sensitively on the details of the heating process, which is 
beyond the scope of this paper, since here we focus on the 
general behavior. In the following, we will assume that the 
temperature is always "jumping" to the highest temper- 
ature solution, which is the solution for minimal cooling, 
and is consequently the one with least heat input. The two 
branches exist over a limited range of redshift, so that this 
ambiguity has no consequences on the modeling. The tem- 
perature behavior is presented in figure ^. As shown there, 
the temperature rises above 10 4 K, between redshifts 5 and 
15, depending on the parameters of the model and reaches 
peak values of the order of a few 10 5 K, that are in agree- 
ment with the Gunn-Peterson test. Interestingly enough, 
in the redshift range 2 to 5, the temperature of the IGM 
lies in a limited range of values above 10 5 K, rather in- 
dependently of the details of the models. The existence 
of a warm IGM is therefore a generic prediction of our 
scenario. 

The instantaneous fraction of cooling gas, / c , is ex- 
pected to be turned into stars after some time and there- 
fore should more or less represent the star formation rate. 
In those models the cooling is spread over a wide range 
of redshifts and starts at rather high redshift. Due to the 
large suppression of the SFR this scenario already pos- 
sesses an interesting property: it leads to an integrated 



amount of stars by the present day which is in reasonably 
good agreement with observations as soon as the feedback 
mechanism is efficient at a level of 10% (i.e. that 10% of 
the bulk of energy produced by supernova is injected in 
the IGM). 

. In such a model the overcooling problem is therefore es- 
sentially solved. This is a considerable success, given the 
limited number of free parameters of the model. However, 
the star formation rate is found to be more or less constant 
with redshift, something which is not in good agreement 
with recent data from the CFRS and HDF, which rather 
indicate late star formation (Madau, 1996). 

4-4- Effects of combined photo-ionization and feedback 

Although we have shown that photo-ionization is clearly 
not a sufficient mechanism by itself to significantly 
limit galaxy formation (at least in order to solve the 
overcooling problem in a realistic scenario), nevertheless 
photo-ionization may play a significant role in the history 
of cooling process (Efstathiou 1992). Furthermore, it 
is clear that photo-ionization is effective from a pure 
observational point of view, since the existence of quasars 
already provides a significant source of UV photons at 
high redshift, even if they might not be entirely respon- 
sible for the ionization. We have therefore examined the 
case of an IGM which would be photo-ionized, and heated 
by some other mechanism (related to star formation) at 
the same time. The effect of photo-ionization is treated 
as in section 3.2, in which we use a simple cooling curve, 
assuming the ionization to be efficient enough to have 
completely suppressed collisional line cooling. 

The temperature has been computed as in the previ- 
ous case. As photo-ionization renders the cooling curve 
regular, essentially by suppressing the line cooling due to 
bound atoms, bifurcations do not appear anymore. The 
resulting temperature history is given in figure |^. The 
behavior of the temperature is more regular than in the 
previous case, essentially due to the regular shape of the 
cooling curve. The temperature at high redshift is similar 
to the pure heating case (previous section), with value of 
the order of a few 10 5 K, insensitive to the details of the 
model. However, there is one noticeable difference at low 
redshift : the temperature is decreasing faster, reaching 
values well below 10 5 K. In our model, the temperature of 
the IGM at some epoch is directly indicative of the mini- 
mum mass of forming galaxies at this epoch. The fact that 
the temperature reach lower values in the present case (see 
figure compared to the pure heating case in figure ||) is 
indicative that in a warm photoionized IGM formation of 
smaller galaxies is allowed if photo-ionization is efficient 
contrary to what one could naively expect. 
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Fig. 7. The integrated amount of stars in term of the 
critical density of the universe is given in two T— CDM 
pictures with T = 0.5 (thin lines) and T = 0.25 (thick 
lines) with different values of the feedback parameter e in 
the photo-ionized case, and fib = 0.05. 
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Fig. 8. An example of the equilibrium temperatures of the 
IGM in the self-regulated picture for the case where the 
IGM is photo-ionized with an additional heating source 
for two specific models : the parameter e in two T— CDM 
pictures with L = 0.5 (thin lines) and T = 0.25 (thick 
lines). The baryonic content of the universe is set to 0.05 
and the normalization of the power spectrum is b = 1.5. 

5. Comparison with observations 

The essence of our model is that galaxy formation is sup- 
pressed at high redshift, because the IGM is pre- heated 
before collapse and that this warm IGM is maintained 
hot by self-regulation with ongoing star formation in 
early galaxies. The most direct observational tests of 
this scenario would be provided by a measure of the 
temperature and the baryonic content of the IGM at high 
redshift. The best measurement of the Gunn-Peterson 
depth at high redshift (Giallongo et al. 1996), suggests 



that a photo-ionized IGM containing most of the baryons 
predicted by nucleosynthesis is ruled out. Numerical 
simulations have however lead to a somewhat different 
picture in which most of the baryons lie in the Lyman-a 
systems. As we have discussed previously, these numerical 
simulations are probably limited by resolution. It is likely 
that dense regions form at high redshift which will still 
exist in the photoionized regime and in the reheated 
area and could then well lead to the Lyman-a systems. 
This scenario has to be elaborated to see whether it can 
actually work, but this is beyond the scope of this paper 
and will be the subject of a future paper. Limits on the 
Compton y-parameter can also provide an interesting 
upper limit on the pressure of the IGM, but are still 
not stringent enough to constrain our scenario. A firm 
prediction of our model is that at redshift ^2 — 4 galaxies 
with circular velocity lower than 100 km/s do just not 
form, although they can form at higher and lower redshift . 



5.1. Neutral gas at high redshift 

The observations of quasars at high redshift have revealed 
the existence of damped Lyman-a systems, which are 
likely to be dense clouds of relatively cool gas. These 
clouds are believed to be the progenitors of present day 
disk galaxies. These clouds can be interpreted in our model 
as transient structures before they commence forming 
stars. Let t g be the characteristic time of survival of the 
gas, if this time scale is short compared to the age of the 
universe tu{z) at epoch z, then the amount of HI gas can 
be related to g(z) via: 



to 



HI 



t H (z) 



g{z)n 



(15) 



In most scenarios, star formation from HI gas is 
triggered by gravitational encounters, this would thus 
mean that the characteristic time scale t* should be of the 
order of the Hubble time. We have therefore estimated 
57hi assuming the ratio . tg , < to be of the order of one. 
The results are presented in figure Bl Our models exhibit 
the correct qualitative behavior for the evolution of the 
HI content. For the CDM picture with T = 0.5, the 
maximum is located at a redshift slightly higher than 
observed, but with L = 0.25 the shape is in satisfactory 
agreement with the data. 

In order to build a consistent picture we select as 
an optimum case only models which reproduce both 
the correct density of stars, and the observed density 
of HI as a function of redshift (Lanzetta et al. 1995; 
Storrie-Lombardi et al. 1996; Natarajan & Pettini 1997). 
Models which lead to a reasonable fl * are not far from 
reproducing correctly the amount of observed HI. This 
is not surprising, given the fact that the observed HI is 
known to be more or less of the order of what is required 
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redshift Z 



Fig. 9. The expected amount of HI against redshift in 
various models, with SI5 = 0.05. Models and lines are the 
same as in figure [| 

to explain present day stars. It is also important to notice 

1 : 



t H (z) 



that this therefore justifies our assumption 
using a smaller value would lead to an overproduction 
of stars. The photo-ionized case with feedback succeeds 
particularly well to match the data for T = 0.25, while 
the case with T = 0.5 fails to fit. This is of interest since 
this offers a constraint on the power spectrum on scales 
much smaller than usually constrained (from clusters or 
CMB data for instance) , still leading to a preferred value 
of T = 0.25. The pure heating case also reproduces the 
observations reasonably, but not in detail: the amount of 
HI gas at low redshift is over-produced by a factor of two, 
while a substantial amount of HI gas appears at high red- 
shift (z > 4) although this might not be regarded as a real 
problem given the absence of information at such redshift. 

The fact that our model reproduces so naturally the 
observed amount of HI is probably its most successful fea- 
ture: in fact, we found that a T — 0.25 photo-ionized warm 
IGM model reproduces automatically the observations as 
soon as the feedback parameter is tuned to reproduce the 
present-day observed amount of stars, a non-trivial result. 

5. 2. Cosmological star formation history 

Given the success of the model thus far, it is tempting to 
go one step further and see whether the model might re- 
produce other features of the cosmic baryon history. One 
of the most important advances in recent years has been 
the estimation of the cosmic star formation history based 
on the CFRS and HDF data. Although star formation is 
certainly a very complex process it is possible that sim- 
ple assumptions may lead to a reasonable picture. Semi- 
analytical models have attempted a first step in this direc- 
tion, but based on a rather different approach than ours. 
A realistic picture would need a more detailed inventory, 



Fig. 10. The amount of HI against redshift in models 
with parameters selected to reproduce quite well the ob- 
served f2» with r = 0.25 in the warm photo-ionized case. 
Dotted, solid and dashed lines correspond respectively to 
the models with (O b ,e) equal to (0.05,0.125), (0.1,0.5) and 
(0.2,2.0). The T = 0.5 case never fit properly the data. 
Data points are taken from Rao & Briggs 1993 (filled cir- 
cle), Natarajan & Pettini 1997 (filled triangles), Lanzctta 
et al. 1995 (filled squares), and Storrie-Lombardi et al. 
1996 (open squares). 
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Fig. 11. Same as figure [h], but without photo- ionization. 
The r = 0.5 case never fits properly the data. 



with attention paid to various types of galaxies, merging 
trees, etc. In order to investigate the star formation his- 
tory in our model, we have therefore adopted a simple rule 
: the fraction of gas available at one epoch is assumed to 
lead to star formation during a fixed period of time At*. 
This time scale represents the typical time scale for con- 
sumption of the gas that will be turned into stars, which 
could be somewhat different from the survival time scale of 
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Fig. 12. Star formation rate estimated by Eq. [l6| for two 

value of the time delay for gas consumption into stars. 
The thin lines are for At* = lGyr, and the thick lines 
are for Ai* = 2Gyr. The data points are taken from Lilly 
et al. 1996 (filled circles), Madau et al. 1996 and 1997 
(filled pentagons), Steidel et al. 1998 (crosses), Conolly et 
al. (open squares), Gallego et al. (filled triangle), Tresse 
& Maddox 1998 (open triangle) and Treyer et al. 1998 
(open circle). The parameters used are those already used 
in figure |l0[ 

the HI gas (t g ) that we introduced previously. The cosmic 
star formation rate at an epoch is therefore: 

P ~ T~F~\ / n b g(z)dt (16) 

We have considered two values of Ai*, 1 and 2 Gyr. 
The resulting cosmic star formation history is plotted in 
figure |lj. Good agreement is found in the low redshift 
regime. This agreement is closely connected to the inclu- 
sion of photo-ionization. The pure heating case does not 
match the data well. Although the modeling of cosmic 
star formation history is likely to be much more complex 
than the simple ansatz we have followed, it is neverthe- 
less interesting to see how well the prediction matches the 
observations. 

5.3. Metal enrichment 

The subject of metal enrichment of the Intra-cluster 
Medium (ICM) has been investigated by numerous 
authors (e.g. see Renzini et al. 1993; Renzini 1997; 
Mushotzky et al. 1996; Ishimaru & Arimoto 1997; Gibson 
et al. 1997; Fukazawa et al. 1998); and it has been advo- 
cated that the ICM could be representative of the IGM as 
far as metals are concerned (Renzini 1998). Self-consistent 
treatments of the metallicity of the IGM in spectrophoto- 
metric models of galaxies are just in their infancy (Sadat 
et al. 1999). As the quantity of iron produced by SNII is 
proportional to the energy released in the IGM, we can 



estimate the abundance of iron from SNII in the IGM in 
our model: 

Zsnii(z = 0) = 0.5 x B*(z = 0)p(M Fe )/(M S Nii) (17) 

where B*(z = 0) is the present fraction of baryons turned 
into stars, (Mp e ) is the average iron mass produced by a 
typical SNII explosion, and the factor 0.5 accounts for an 
equipartition of the iron between the stars and the IGM. 
Adding the contribution of the SNIas which provide 75% 
of the iron in the IGM (see Renzini et al., 1993):we get: 

Z t ot(z = 0) = Z SNII + Z SNIa = 4 x Z SNIIs (z = 0) (18) 

One thus finds typical metallicities of the IGM of 0.1 so- 
lar at redshift z = 0, and of 0.05 solar at redshift z = 2.5, 
which are compatible with recent observations in Damped 
Lyman a Systems (Lu et al. 1996; 1997), these systems are 
assumed to be representative of the state of the IGM as 
far as metals are concerned. This is not a strong constraint 
on our model however, as the metalicity of the IGM essen- 
tially reflects the metal production inferred from the star 
formation rate. 

5.4- Discussion and Conclusion 

We have examined in detail the overcooling problem once 
again. We confirm the previous claim that in the absence 
of significant feed-back mechanism, cold gas, eventually 
in the form of stars, is over-produced in Cold Dark Mat- 
ter models. Observations reveal a little amount of stars 
compared to standard nucleosynthesis prediction. A fun- 
damental question is therefore: where are the primordial 
baryons ? There are reasonable arguments which indicate 
that they are in the intergalactic gas, maintained at some 
temperature above 10 4 K. We found that photo-ionization 
reduces significantly the amount of cooling, especially in 
the case of a low baryonic density. However, from our anal- 
ysis, we found that the photo-heating is not sufficient to 
suppress the cooling of the gas, over-producing in a signif- 
icant way the observed amount of cold gas. This strongly 
argues for the existence of strong feed-back mechanisms 
during galaxy formation. We have explored in some de- 
tail the warm IGM picture in which the high redshift gas 
is heated preventing further galaxy formation. We found 
that the IGM can be easily heated to temperatures of 
the order of 10 5 — 10 6 K, altering substantially the pro- 
cess of galaxy formation. In such models, the overcool- 
ing problem is easily solved provided that at least 10% 
of the energy of supernova is transferred to the IGM. We 
found that with this hypothesis the observed amount of 
stars can be reproduced for very reasonable values of the 
feedback parameter. Such scenarios reproduce quite well 
the observed amount of HI gas at high redshift, especially 
when photo-ionization is taken into account. This is a re- 
markable success of this model. Assuming that this gas is 
turned into stars, we also found that the cosmic stars for- 
mation rate is well reproduced. Clearly, this model seems 
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to meet success, despite the fact that the galaxy formation 
history is significantly different than in more traditional 
scenario in which large galaxies are built from the merg- 
ing of smaller entities. This certainly illustrates the fact 
that large galaxies may well have been formed at relatively 
high rcdshift (say between 3 and 5), in T-CDM model with 
r = 0.25. High rcdshift galaxy formation of large galaxies 
is not therefore intimately linked to low density universes 
(although very high redshift - greater than 10 - would cer- 
tainly be a problem in a high universe). 

It is interesting to notice at that point that this 
scenario provides us also with a natural explanation for 
the Lx — Tx departure from scaling law for small clusters 
of galaxies (Kaiser 1991; Evrard & Henry 1991). Indeed, 
if the IGM is pre-heated before the formation of the 
clusters, then the zero-point entropy of the gas prevents 
it to be very concentrated in the cores of clusters. This 
effect is very important when this zero-point entropy is 
comparable or bigger than the entropy created by shocks 
during the formation of the cluster, and then affects pri- 
marily small galaxies, clusters and groups. This departure 
from pure scaling can also be seen in the evolution of the 
surface brightness profiles for different gas temperature 
ranges (data taken from ROSAT and GINGA, Ponman et 
at 1998). These authors estimate this zero-point entropy 
to be of the order of 100 /i -1 / 3 keVcm 2 , which can be 
provided by a uniform IGM pre-heated to a temperature 
Tjgm ~ 310 4 (1 + z) 2 K. The required temperature is 
then ~ 2.5 10 5 K at z = 2, which is compatible with the 
results of our model (see Fig. |J). 

In our scenario, galaxy formation occurs in two clearly 
distinct phases. At high redshift objects which form are 
cooling with a cooling time much shorter than the age of 
the universe. This could be argued as being the phase of 
bulges and elliptical formation. In this regime, the forma- 
tion history follows the classical hierarchical picture, small 
objects form first, large objects latter. At later epochs, 
this scheme is reverted. This might well be the epoch of 
disk formation, during which the gas gently falls in some 
potentials which already contain stars from the previous 
generation. This would imply that the typical epoch for 
disk formation is z = for the smallest to z — 2 — 3 for 
the largest disks. Beyond this epoch most of the form- 
ing structures would correspond to bulges and ellipticals. 
This scheme is certainly consistent with observations as 
high redshift objects seem to have characteristics of bulges 
in their early stages. A firm prediction of our model is 
that galaxy formation with circular velocity smaller than 
100 km/s should be strongly suppressed at redshift in the 
range 1 — 5. Clearly, observations of the dynamical state 
of high redshift galaxies would be important to test this 
scenario. This is probably testable in future NGST ob- 
servations. Of course, the most direct test of our scenario 
would be to measure the temperature history of the IGM. 
However, this seems to be difficult. 
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